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Nonleptonic two body Be decays including radially excited 'tp{2S) or rjdfiS) mesons in the final 
state are studied using the perturbative QCD approach based on kx factorization. The charmonium 
distribution amplitudes are extracted from the n = 2,1 = 0 Schrddinger states for the harmonic 
oscillator potential. Utilizing these distribution amplitudes, we calculate the numerical results of 
the Be —^ 'if{2S),rie{2S) transition form factors and branching fractions of Be —^ 'if{2S)TT,'ge{2S)n 
decays. The ratio between two decay modes Be —>■ iIj{2S)h and Be —>■ J/inT is compatible with the 
experimental data within uncertainties, which indicates that the harmonic-oscillator wave functions 
for 4’{‘2S) and r]e[2S) work well. It is found that the branching fraction of Be rie{2S)n, which 
is dominated by the twist-3 charmonium distribution amplitude, can reach the order of 10~®. We 
hope it can be measured soon in the LHCb experiment. 

PACS numbers: 13.25.Hw, 12.38.Bx, 14.40.Nd 


I. INTRODUCTION 

The meson B^, a pseudoscalar ground state of h and c quarks, can only decay through weak interactions. Either 
of the heavy quarks (6 or c) in it can decay individually, which makes it an ideal system to study weak decays of 
heavy quarks. Around 0(10®) mesons can be anticipated with 1 fb“^ of data at the LHC [I|, which is sufficient for 
studying the meson family systematically. Up to now, several new decay channels of the B^ meson have 
been successfully observed by the LHCb Collaboration, while an excited Be meson state which is consistent with 
expectations of the Bc{2S) has been found by the ATLAS detector Q. 

Recently, the LHCb Collaboration observed the decay mode Be —>■ iP{2S)tt for the first time with the measured 
ratio of the branching fractions as Q 

= 0.250 ± 0.068(stat) ± 0.014(syst) ± 0.006(H) . (1) 

B{Bc J/lpTT) 

The last term above accounts for the uncertainty on B{ijj{2S) —?> yL~)/B{J/ip —)• pL~). Although there is not much 

data for the Be meson decaying into two-body final states containing a radially excited charmonium such as ip(2S) 
or ijc{2S) except the Be —>■ ip{2S)Tr channel, many theoretical studies of nonleptonic Be decays with radially excited 
charmonium mesons in the final state have been performed by using various approaches. For example, in Ref. Q, 
the authors computed the branching ratios for Be —^ 'tl){2S)X decays with the modified harmonic-oscillator wave 
function in the light front quark model; in Ref. [l^, the ISGW2 quark model was adopted to study the production of 
radially excited charmonium mesons in two-body nonleptonic Be decays, the relativistic (constituent) quark model, the 
potential model, the QCD relativistic potential model , an d the improved instantaneous BS equation and Mandelstam 
approach were adopted in Refs, [ill [Hi, R-ef- [13 > 13 and Ref [3j respectively. However, all of these calculations 

are based on a naive factorization hypothesis, with various form factor inputs. There are uncontrolled large theoretical 
errors with quite different numerical results, and most of them cannot give any theoretical error estimates because of 
the unreliability of these models. 

The perturbative QCD approach (pQCD) [3 based on kp factorization, which not only can deal with the emis¬ 
sion diagrams corresponding to the naive factorization terms basically, but can also handle well the nonfactorizable 
diagrams by introducing the wave function of the light meson in the final states of the Be decay modes, is widely used 
in the nonleptonic two-body Be decays [TtI - I^ • In onr recent work [3 , the pQCD approach was used successfully in 
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describing the S-wave ground state charmonium decays of B,. meson based on the harmonic-oscillator wave functions 
for the charmonium IS states. In this work, we will use the harmonic-oscillator wave function as the approximate 
wave function of the 2S states and study the Be —>■ i/^( 2 S') 7 r, r]c{2S)n decays in the pQCD approach to provide a ready 
reference to the existing and forthcoming experiments. 

The structure of this paper is organized as follows. After this Introduction, we describe the wave functions of 
radially excited charmonium mesons tf}{2S),r]c{2S) in Sect. II. We calculate and present the expressions for the 
Be —>■ tp{2S), f]e{2S) transition form factors in the large-recoil regions and the Be —>■ ^/>( 2 S') 7 r, r]e{2S)n decay amplitudes 
in Sect. III. The numerical results and relevant discussions are given in Sect. IV, and Sect. V contains a brief summary. 


II. WAVE FUNCTIONS 

In hadronic B decays, there are several energy scales involved. In the expansion of the inverse power of heavy quark 
mass, the hadronic matrix element can be factorized into perturbative and nonperturbative factors. In the pQCD 
approach, the decay amplitude A{Be M 2 M 3 ) can be written conceptually as the convolution [I^ 

AiBe ^ M 2 M 3 ) - Jd*kid'^k 2 d*k 3 TT[C{t)<^>BAkl)^M,ik2)^M,ik3)H{ki,k2,k3A)] , (2) 

where fc^’s are momenta of spectator quarks included in each meson, and “Tr” denotes the trace over Dirac and color 
indices. In the above convolution, C{t) is the Wilson coefficient evaluated at scale t, the function H{ki,k 2 ,k 3 A) 
describes the four quark operator and the spectator quark connected by a hard gluon, which can be perturbatively 
calculated including all possible Feynman diagrams without end-point singularity. The wave functions $m 2 

and $m 3 describe the hadronization of the quark and anti-quark in the Be meson, the charmonium meson il’{2S) or 
rie{2S) and the hnal state light meson pion, respectively. 

As a heavy quarkonium discussed in Refs. |^, the nonrelativistic QCD framework can be applied for the Be 
meson, which means its leading-order wave function should be just the zero-point wave function with the distribution 
amplitude 

~ ^c/mBjexp[-uj%y/2] . (3) 

For the light meson pion, we adopt the same distribution amplitudes and as defined in Refs. [sTI.!^. 

The harmonic-oscillator wave functions has been adopted to describe the IS state mesons [33l - [^ . and they can 
explain the experimental data well [2^. In the quark model, r]e{2S) and ip{2S) are the first excited states of r]e and 
J/')/', respectively. The 2S means that for these states, the principal quantum number n = 2 and the orbital angular 
momentum I = 0. The definitions of the 2S state wave functions are similar to the IS states via the nonlocal matrix 
elements [s^, 

(V'( 25 ')(P,e'^)|c(z)„c( 0 ); 3 | 0 ) = ^ 

{^{2S){P,e^)\c{z)e,c{0)p\0) = + {f 

{r]e{2S){P)\c{z)ac(0)/s\0) = --^^= j da:e“-^'*[( 75 /^)a/ 3 ' 0 ”(a;, 6 )-hm( 75 )„; 3 ' 0 ®(x,&)], (4) 

where P stands for the momentum of the charmonium meson rje(2S) or ip(2S) and m is its mass. The x represents 
the momentum fraction of the charm quark inside the charmonium, and b is the conjugate variable of the transverse 
momentum of the valence quark of the meson. The AA) denotes its longitudinal (transverse) polarization vector. 
The asymptotic models for the twist-2 distribution amplitudes and the twist-3 distribution amplitudes 

will be derived following the prescription in [s^. 

First, we write down the Schrddinger equal-time wave function for the harmonic-oscillator potential. The 

radial wave function of the corresponding Schrddinger state is given by 

^'( 2 S)(r) cx (^ - (5) 

where A = and w is the frequency of oscillations or the quantum of energy. We perform the Fourier transformation 
to the momentum space to get 'I' 2 s(k) as 



'I'( 2 S)(k) oc (2fc^ - 30^)6 


( 6 ) 
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with being the square of the three momentum. In terms of the substitution assumption, 


ku^kj_, kz^^{x-x), ml = ^ , 


XX 


with rUc the c-quark mass and x = 1 — x. We should make the following replacement as regards the variable 

,2 . k]_ + {x-xfml 

* ^ • 


( 7 ) 

( 8 ) 


Then the wave function can be taken as 


^( 2 S)(k) ^'( 2 S)(a;,k_L) oc ( 


k5_ + ml(x — x)^ 
2xx 




— 3o; )e Sxxa 


(9) 


Applying the Fourier transform to replace the transverse momentum kj^ with its conjugate variable b, the 2S oscillator 
wave function can be taken as 


^(2S)(a::,b) - J (fk±e ku) 




with 


T{x) = 1 — Ah^mcUJXx + 


mc{x — xY 


OJXX 


We then propose the 2S states distribution amplitudes inferred from Eq. (US, 


( 10 ) 


( 11 ) 


( 12 ) 


with the $“'*^(a;) being the asymptotic models, which have been given in [s^- Therefore, we have the distribution 
amplitudes for the radially excited charmonium mesons r]c(2S) and '0(25') 


,L,T,V/^ _ hs 




2C27Ve 


vT{x)e 




'^Yx,b) = 


2C27Ve 


N*{x — xYT{x)e 


«'^(a;, 6 ) = 

^'"(a;, 6 ) = 

with the normalization conditions: 


V(^ M _ hs atV!^ , + 


2C27Ve 

hs 

2y/Wl 


[1 + (a; — xY]T{x)e 


N^T{x)e 




(13) 



FIG. 1: The shape of the distribution amplitude for ip^{x) when 6 = 0, with the solid (dashed) line for ix = 0.2(0.3) GeV. 
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4>*(x, 0)dx 
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(14) 


Nc above is the color number, N'^(i = L,T^t,V,v, s) are the normalization constants, and / 2 s is the decay constant 
of the 2S state. All the distribution amplitudes in Eq. (TT^ are symmetric under x O Here we do not distinguish 
the leading twist distribution amplitude of the r]c{2S) meson from 4/^^^ of the '0(25') meson, and the same decay 
constant has been assumed for the longitudinally and transversely polarized 0(25) meson. To make things clearer, 
the shape of the distribution amplitude 'I'^(x, 0) is displayed in Fig. [T] The free parameter uj = 0.2 GeV is adopted, 
such that the valence charm quark, carrying the invariant mass ~ m^, is almost on shell. It can be seen that 

the two maximum positions are near x = 0.35 and x = 0.65 and a larger value of parameter w gives a wider shape. 
Note that the dip at x = 0.5 is a consequence of the radial Schrddinger wave function of the n = 2,1 = 0 state. 


III. FORM FACTORS AND DECAY AMPLITUDES 
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FIG. 2: The leading-order Feynman diagrams for the Be —>■ (0(25'), 77 c(25')) transitions. 


In the pQCD approach, the Be —t tjj(2S),r]c{2S) transition form factors in the large-recoil limit (g^ = 0), which are 
similar to that of Be —>■ J/tp,rie [1^, can be calculated from above universal hadronic distribution amplitudes. The 
lowest-order diagrams are displayed in Fig. [5J The form factors V{q^) and Aq, 1 , 2 ( 9 ^) are defined via the 

matrix element [39|, 


{rie{2S)[P2)\crb\Be{Pi)) 

(0(25)(P2)|c7'^&|Sc(Hi)) 

mS){P2)\criBh\Be{Pl)) 


{Pi + P2Y - 




F+{qY + 




q^Fo{qY, 




e:P2pPla , 


M + m 

2mAo{qY^^q^ + {M + m)Ai(g^) 




M + m 


{Pi + P 2 Y - 


(15) 

(16) 

(17) 


where q = Pi — P 2 is the momentum transfer and PYP 2 ) is the momentum of the initial (final) state meson. M is the 
mass of Be meson, and e* is the polarization vector of the 0(25) meson. In the large-recoil limit, say q^ = 0, we have 

Fo(0)=F+(0), Ao(0) = i^Ai( 0 )-i^A 2 ( 0 ). (18) 


It is straightforward to calculate the form factors Fo(g^), V{qY and Ao,i(g^) at the tree level in the pQCD. They read 
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i^o = 2 <J^TrM'^fBCf j dx2 j 6162^61^62 exp 


2 

X I [i/;’'(a; 2 , 62 ) (x 2 - 2n) - 'ip%X 2 ,b 2 )r (2x2 - n)] Eab(ta)h(ae, /3a, h,b2)Stix2) 

+ [■;/'”(x2,62)(rc + r^(l - xi)) - i/>®(x2, 62)2r(l - Xi + Vc)] Eab{tb)h(ae, l 3 b,b 2 ,bi)St{xi) L 


( 19 ) 


y = 2 y-(l + r)TrM'^fBCf J dx2 J 6162^611/62 exp (-yJ-j 

X I [^^(X 2 , 62 )r (1 - X 2 ) +'tlj'^{x2,b2)(rb - 2)] Eab(ta)h(ae, l3a,bi,b2)St(x2) 

- {x 2 ,b 2 )rEab{tb)h(ae, /3b, 62, 6i)S't(xi) 


( 20 ) 


^0 = 2 Y^ 7 rM^/BC'/ j dx 2 J 6i62d6id62exp 

X I [■i/’'^(x2,62) (X2 - 2 rb) - ^*(x2,b2)r (2X2 - rb)] Eab{ta)h(ae, /3a, 61,62)S't(X2) 
- ijj^(x 2 ,b 2 ) [tc + r^(l - Xi)] Eab(tb)h(ae, l3b,b2,bi)St(xi) i. 


( 21 ) 


Ai = 2 


2 r 


31 + r 


^CXD 

ttM^/bC/ / (ix 2 / bib2dbidb2 exp 

Jo Jo 


ulbl 


'ip^{x2, 62) (1 + X2 - r^(l - X2) - 4 rb) + V'^(x2,62) 


r(2 - 4x2 + n) + 


Tb - 2 


X Eab{ta)h(ae, l3a,bi,b2)Stix2) - (2^2,62) [l “ 2 xi + 2rc + Eabitb)h{ae, Pb,b 2 ,bi)Stixi) L ( 22 ) 


with r = ^ and rb,c = The functions Eab, the scales ta,b and the hard functions h are given in Appendix B of 

Ref. m. 

The quark diagrams contributing to the Be —>■ iIj{2S)tt, r]c{2S)TT decays are displayed in Fig. [31 where (a) and (b) are 
for the factorizable topology, and (c) and (d) are for the nonfactorizable topology. The effective Hamiltonian relevant 
to the considered decays is written as [i^ 

-Heff = ^V:,Vu4Cli^i)0^(^,) + C2i^i)02(^^)] + h.c., (23) 

with V*ij and Vud the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, Ci_ 2 (/r) the Wilson coefficients, and 
Oi, 2 (At) the effective four quark operators 

Ol(/i) = 6a7''(l - 75)C/3 0 M“/37ij( 1 - 75)da, 

02(/i) = 6a7'"(l - 75)Ca O W/37 m( 1 - 75)d/3, (24) 

where a and /3 are the color indices. Since the four quarks in the operators are different from each other, there is 
no penguin contribution. Therefore there will be no CP violation in the decays of Be —> '0(2S')7r, ?7c(2S')7r within 
the standard model. After a straightforward calculation using the pQCD formalism of Eq. ([3]), we have the decay 
amplitudes 

AiBe ^ {4(2S),r)e(2S))7r) = V4Vud[(C2 + \ci)Te + CiMe]- (25) 

The detailed expressions of Te and A4e are the same as the Be —>■ (J/i/’, Vc)'^ decay modes in Appendix A of Ref. [^ . 
except for the replacements J/ijj ^ 'ij){2S) and rje ^ r]e{2S). 
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FIG. 3: Feynman diagrams for Be —>■ 'tp{2S)n,rjc{2S)Ti decays. 

IV. NUMERICAL RESULTS AND DISCUSSIONS 


In the numerical calculations we need the following input parameters (in units of GeV) [4l|: 

rric = 1.275, mb = 4.18, Mb^ = 6.277, m^( 2 S) = 3.686, = 3.639. (26) 

For the relevant CKM matrix elements we use Vcb = (40.9 ± 1.1) x 10“^ and Vud = 0.97425 ± 0.00022 [4ll| . 

The decay constant f 11 ,( 23 ) can be derived from the process e+e“ by the relationship 


using the data given in 


im 


/)/>(2S) — 


I ^mii,(2S)^il,(2S)—ye+e 




r,/,( 2 S)-)-e+e- = (2.36 ±0.04) keV. 


(27) 


(28) 


Then we have fip( 2 S) = 296 I 2 hleV. The decay constant fne( 2 S) can be determined by the double photon decay of 
r]c{ 2 S) as 


/r,c(2S) - 


I 8l7r?7i^^^2S)4^?7c(2S)—>77 


y 4(47ra)2 

Using the measured results of the branching fractions ric{2S) —?> 77 and the full width of r]c{2S) [4ll| . 


(29) 


Bir)c{2S) ^ 77 ) = (1.9 ± 1.3) x 10-^ r^^( 2 S) = 11 , (30) 

we can get the decay constant fne( 2 S) = 2431™^ MeV. As for the decay constant for Be, we adopt /b_, = 489 MeV (i^ . 

Our numerical results for the form factors yBc^ip( 2 S) listed in Table HI We find 

that the form factors are close by different approaches within errors, except the results in Ref. which are typically 
smaller. Some dominant nncertainties are considered in our numerical values: the first error comes from the shape 
parameters ujb = 0.6 ± 0.1 (w = 0.2 ± 0.1) GeV for the Bc{'4>{2S)/r]c{2S)) meson, the second one is induced by 
me = 1.275 ± 0.025 GeV, the third error comes from the decay constants of the ip{2S) or ric{2S) meson, and the last 
one is caused by the variation of the hard scale from 0.75t to 1.25t in Eq. which characterizes the size of next- 
to-leading-order contribution. It is found that the main errors come from the uncertainties of the shape parameters 
and the charm-quark mass. Therefore, the decay of Be —>■ ip{2S){rie{2S)) provides a good platform to understand the 
wave function of the radially excited charmonium states and the constituent quark model. The uncertainty from the 
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TABLE I: The form factors for ^b^( 2 S) yB,^i,( 2 S) ^ 0 evaluated by pQCD and by other methods 

in the literature. We also show theoretical uncertainties induced by the shape parameters, rUc, /t/;( 2 s) or fr]c( 2 s) ^'iid the hard 
scale t, respectively. 



This work 

Ref. 19] 

Ref. 110]“ 

Ref. Ill] 

Ref. 143] 

To 

Ao 

Ai 

A 2 

Vo 

^ 7r.+0.09+0.12+0.23+0.02 
U- <U_o.05-0.10-0.32-0.01 
^ r^+0.09+0.07+0.00+0.01 
U.OD_o.05-0.04-0.00-0.01 
n c^+o.13+0.06+0.00+0.01 
U.OO_o.04-0.03-0.00-0.01 
/A f^9+0.27+0.04+0.01+0.02 
^•^^-0.05-0.01-0.01-0.00 
f. Qc+0.18+0.15+0.01+0.03 

u.yo o. 08 - 0 . 10 - 0 . 01 - 0.01 

0.45 

0.335 

0.102 

0.525 

0.325 

0.42 

0.35 

0.15 

0.73 

0.27 

0.23 

0.18 

0.14 

0.24 

0.20 

0.38 

0.90 

0.90 


“Comparing the definitions of the transition form factor of Ref. [Toll with ours, we have the following relations at the maximal recoil 
point: 

To = /+, V = {M + m)g, = —^ -, A 2 =-(M + m)a+, 

M + m 

f + (M^ — m^)a^ + q'^a- 

Ao = ---, 

2m 

where the values of /+ ,g,f, a+, a— are given in 0- 


TABLE II: Branching ratios (10 of the Be -A r)c{2S)'K,'ip(2S)'R decays. The errors induced by the same sources as they in 
Tabled! 


Modes 

This work 

[9]“ 

110] [11] 

112] 

113] 


[15] 

Be —>■ rje(2S)'!T 
Be — >■ 'tp{2S)'K 

-.r, O+3.4+4.0+7.8+1.2 
-L0>.O-i_8-2.8-7.2-0.4 
^ 7+2.8+1.8+0.1+0.7 
'-1.1-1.2-0.1-0.3 

2.97 

2.4 1.7 

3.7 1.1 

2.2 

0.63 

2.4 

2.2 

0.66 

2.0 

2.87 

2.66 7.6(5.8) 


“We quote the result with the modified wave functions for -(/)(25). 

^The nonbracketed (bracketed) results are evaluated at the NLO (LO) level. 


decay constant of r]c{2S) meson is large due to the low accuracy measurement of the branching fraction in Eq. (1301) : 
the relevant uncertainty of Fq is large, too. We expect that it could be measured precisely at LHCb and Super-B 
factories in the near future. We also noticed that the error from the uncertainty of the hard scale t is small, which 
means the next-to-leading-order contributions can be safely neglected. The errors from the uncertainty of the CKM 
matrix elements are very small and they have been neglected. 

The branching fractions for the Be -A ric(2S)Tr,tjj(2S)TT decays in the meson rest frame can be written as 


S(F,^ (^(2F),ryc(2F)K) = 


32^ Mb 




(31) 


where the decay amplitudes A have been given explicitly in Eq. (1251) . In Table HIl we show the results of the branching 
fractions for the two-body nonleptonic Be —f rie(2S)n, 'iI)(2S)tt decays, where the sources of the errors in the numerical 
estimates have the same origin as in the discussion of the form factors in Table H] It is easy to see that the most 
important theoretical uncertainties are caused by the nonperturbative shape parameters, the charm-quark mass, and 
the decay constant /,,„( 2 S), which can be improved by future experiments. It is found that the branching fractions 
of Be decays to the 2S state are smaller than those of IS state in our previous study [ 2 ^ in the perturbative QCD 
approach. This phenomenon can be understood from the wave functions of the two states. The presence of the node in 
the 2S wave function, which can be seen in Fig.[Tl causes the overlap between the initial and final state wave functions 
to becomes smaller. Besides, the tighter phase space and the smaller decay constants of 2S state also suppress their 
branching ratios. 

We also make a comparison of our results with the previous studies. One can see that our results are comparable to 
those of [ 4 ^ within the error bars, but larger than the results from other modes. This is because they have used the 
smaller form factors at maximum recoil. Regardless of this effect, our results are consistent with theirs. For example, 
as shown in Tables U and im our values of Aq and Fq are about 2.5 times the results of Ref. [ll|, and result in our 


branching ratios are 6 times larger than theirs. For a more direct comparison with the available experimental data, 
we compare the present results in Table m with those for the decays of Be to S-wave charmonium states J/ij} and rje 
(also based on the harmonic-oscillator wave functions), whose results can be found in Ref. [ 2 ^ , and obtain the ratios 
B{Be —t {'ijj{2S)TT))/B{Be —t {J/ijjTr)) = 0.291q and B{Be —t (r]e(2S)Tr))/B(Be —t {rje'K)) = O.SSIqJq. The former is 





















TABLE III: The values of decay amplitude from twist-2 and twist-3 charmonium wave functions for Be —>■ rie{2S)'n',ip{2S)n 
decays. The results are given in units of GeV®. 


Modes 

twist-2 

twist-3 

total 

A{Bc -t iji2S)n) 

-1.7-0.07i 

-0.4-0.06i 

-2.l-0.13i 

A{Bc —t r]c{2S)Tv) 

-1.5-2.3i 

3.9-bl.4i 

2.4+0.9i 


consistent with the data 0.25 ±0.068 ±0.014 Q, and also comparable with the recent prediction of the Bethe-Salpeter 
relativistic quark model [l^, 0.24. This fact may indicate that the harmonic-oscillator wave functions for radially 
excited states are reasonable and applicable. Although the Be —>■ ?7c(2S')7r decay has not yet been measured so far, 
the predicted large branching ratio (10“^) makes it possible to measure it soon at the LHCb experiment or future 
facility. 

We now investigate the relative importance of the twist-2 and twist-3 contributions in Eq. (jT]) to the decay amplitude, 
whose results are displayed separately in Table Hill where the label “twist-2 (twist-3)” corresponds to the contribution 
of the twist-2 (twist-3) distribution amplitude only, while the label “total” corresponds to both of the contributions. 
It is found that the contribution of twist-3 distribution amplitude is not power-suppressed for Be —>■ r]e{2S)'!r decay, 
whose contribution is 1.5 times larger than the twist-2 contribution. The reason is that the term (x2,b2)2r in 
Eq. (froi) from Eig. [3] b gives the dominant contribution to the decay amplitude, since the asymptotic model of the 
twist-3 distribution amplitude in Eq. m for the rie{2S) meson has no factor like x{l — x) to suppress its integral value 
in the end-point region, which leads to large enhancement compared with twist-2 contribution. However, because 
twist-3 terms of the ip{2S) meson distribution amplitude do not contribute to the Be —^ 'iP{2S)tt decay amplitude from 
FigElb, the contribution from other diagrams with twist-3 distribution amplitude is only one-fifth smaller than that of 
the twist-2 contribution in this process. It is also found that there is very strong interference between contributions of 
the twist-2 and twist-3 wave functions for both Be —>■ '0(25')7r and Be r]c{2S)'K decays. The numerical results show 
that the contributions from twist-3 wave function have an opposite sign between the two channels. This results in 
constructive interference for the former, but destructive interference for the latter. The reason is that the amplitudes 
are different between the two decays at twist-3 level, which can be seen in Eqs. (Al) and (A4) of [2^. A similar 
situation also exists in Be —>■ D-k, D*tt (2^ decays. 


V. CONCLUSION 

We calculated the form factors of the weak Be decays to radially excited charmonia and the branching ratios of 
Be —>■ ip{2S)Tr,r]e{2S)Tr decays in the pQCD approach. The new charmonium distribution amplitudes based on the 
radial Schrddinger wave function of the n = 2, / = 0 state for the harmonic-oscillator potential are employed. We 
discussed theoretical uncertainties arising from the nonperturbative shape parameters, the charm-quark mass, the 
decay constants, and the scale dependence. It is found that the main uncertainties of the processes concerned come 
from the shape parameters and the charm-quark mass. The theoretically evaluated ratio B{Be —>■ {'tp{2S)TT))/B{Be —>■ 
{J/Pjtt)) = 0.29t°:ll is consistent with the data, which indicates the harmonic-oscillator wave functions work well, not 
only for the ground state charmonium, but also for the radially excited charmonia. It is also found that the twist-3 
charmonium distribution amplitude gives a large contribution, especially for Be r]e{2S)Tr decay, whose branching 
fraction is of the order of 10“^, which could be tested at the ongoing large hadron collider. 
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